Two Dimensional Poynting Flux Dominated Flow onto a Schwarzschild Black Hole 
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We discuss the dynamics of the accretion flow onto a black hole driven by Poynting flux in a 
simplified model of a two-dimensional accretion disk. In this simplified model, the condition of 
the stationary accretion flow is found to impose a nontrivial constraint on the magnetic field con- 
figuration. The effect of the magnetic field on the accretion flow is discussed in detail using the 
paraboloidal and hyperboloidal type configuration for the poloidal structure suggested by Bland- 
ford in 1976. It is demonstrated explicitly that the angular velocity of the disk, Qd, deviates from 
the Keplerian angular velocity. The angular velocity of the rigidly-rotating magnetic surface, SIf, 
does not have to be the same as the angular velocity of the disk for the paraboloidal type configu- 
ration. But for the hyperboloidal type configuration, it is found that we can set Hp = i^D, which 
corresponds to an accretion disk of perfect conductor. We discuss the numerical solutions of the 
stream equation for stationary accretion flow in the Schwarzschild background using a paraboloidal 
type configuration. The dynamics of the accretion disk is found to depend strongly on the ratio of 
the accretion rate to the magnetic field strength. 

PACS numbers; 97.10.Gz, 97.60.Lf 
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I. INTRODUCTION 

To describe the powerful and highly coUimated astro- 
physical jets observed in AGN and quasars, Poynting 
flux model has been suggested long time ago 0, Q and 
many interesting works have been developed. One of the 
characteristics of Poynting flux is that it carries very lit- 
tle baryonic component compared to the hydrodynamic 
flow. Recently this property of Poynting flux is found 
to be consistent with the required property for powering 
GRBJI . And the Poynting flux in a system of black hole- 
accretion disk has also been studied in connection with 
the gamma ray bursts 0, 0, and it is found that the 
evolution of the system is largely depend on the Poynt- 
ing outflow from the disk0, U • 

The configuration of the ordered magnetic field around 
the accretion disk, which is responsible for the Poynting 
flux, has been discussed both in analytical and numeri- 
cal studies. For the strong enough electromagnetic held 
around the compact object, the force-free magnetosphere 
can be established. In the non-relativistic formulation, 
Blandford'^ suggested an axisymmetric and stationary 
solution for the Poynting outflow, assuming a force- 
free magnetosphere surrounding an accretion disk. The 
poloidal field configuration for a black hole in a force-free 
magnetosphere has been discussed recently by GhoshQ 
in the relativistic formulation using Grad-Shfranov equa- 
tion. The developments of the ordered magnetic field in 
the disk and the Poynting outflow from the disk have 
been studied numerically by many authors^, ITll H^ . 
Recently Ustyugova et al.jlSi] performed an axisymmet- 
ric magneto-hydrodynamical simulation to show that the 
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quasi-stationary and approximately force-free Poynting 
jet from the inner part of the accretion disk is possible. 

In studying the accretion flow, it is also important to 
know the dynamics of the accretion disk influenced by the 
magnetic field. The energy and angular momentum car- 
ried out by the Poynting flux are determined essentially 
by the torque and the angular velocity of the magnetic 
surface D, p , which depends strongly on the disk dynamics 
subjected to the magnetic field. We adopt a simplified 
two-dimensional accretion disk model^3| embedded in a 
stationary, axisymmetric and force-free magnetosphere. 
To see the magnetic effect transparently it is also as- 
sumed that there is no viscous stress tensor in the disk 
and there is no radiative transfer from the disk in this 
model. 

The state of an accretion flow is governed by a set of 
equations obtained from the conservation of stress-energy 
tensor 14] and the force-free magnetosphere with the cur- 
rent configuration is determined by the Grad-Shafranov 
equation. The two sets of equations are not completely 
independent to each other because they are coupled at 
the disk on the equatorial plane. Hence the variables in 
the Grad-Shfranov equations, the stream function '5, the 
current potential / and p are constrained by the accre- 
tion equations. In other words, the accretion equations 
provide a kind of boundary values of 5*, / and Qp which 
cannot be determined by Grad-Shafranov equation alone. 
It also impose a nontrivial constraint on the configuration 
type of the magnetic field. For example, it is found that 
the hyperboloidal configuration suggested by Blandford 
in 19760 is consistent with an accretion disk of perfect 
conductor, that is 51^;' = Qpi, for the stationary accretion 
flow, whereas for the paraboloidal type conflguration Q p 
does not have to be the same as fin . 

In section 2, we describe the simplifled two-dimensional 
accretion disk model dominated by Poynting flux in a 
force- free magnetosphere. The system is governed by four 
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basic equations: three accretion equations and the Grad- 
Shafranov equation in the Schwarzschild background. In 
section 3, numerical analysis of the accretion flow and 
the field configuration of paraboloidal type are discussed 
in detail for two cases of lower and higher accretion rates. 
The results are summarized and discussed in section 4. 



II. TWO DIMENSIONAL ACCRETION DISK IN 
FORCE-FREE MAGNETOSPHERE 

In this work, a black hole-accretion disk system is sup- 
posed to be surrounded by the ordered electromagnetic 
field, which is determined by the electric charge and cur- 
rent distributions as a solution of the inhomogcneous 
Maxwell equation: 



F'!,': = 47r 



(1) 



where are the fi-th component of the bulk electro- 
magnetic current density. The boundary values for the 
field configuration are responsible to the surface current 
densities on the horizon and/or on the accretion disk. 
It is a kind of boundary value problem. which can 
be described by introducing appropriate surface current 
densities, j^, defined on the surface such that the total 
conserved current density, is given by 



Using the current conservation, 



(2) 



(3) 



one can identify the surface current density when pro- 
vided with the geometry of the surface. 

The horizon of a black hole is a mathematically well 
defined surface, on which the appropriate surface current 
density has been studied in depth 16]. On the other hand, 
the surface of the physical accretion disk is not expected 
to have any sharp boundary and various shapes have been 
suggested depending on the properties of accretion flow. 
In this work, however, we assume a simplified accretion 
disk, which is vertically squeezed down to the equatorial 
plane, a two-dimensional accretion disk. Then we can 
associate the surface current density in a simple way. 



3^ ^^{Fl^ -F'j')5{e-i,/2), 



(4) 



which are nothing but Gauss' law and Ampere's law as 
given by 



(5) 



= - — (Bi-Bt), = —(Bi- BL), (6) 

An 47r 



where and are surface charge and surface cur- 
rent density(spatial) respectively and +/— denote the 
upper/lower hemisphere. 



A. Accretion equations 

The dynamics of an accretion flow is determined by 
the conservation equation of the stress-energy tensor: 



(7) 



The stress-energy tensor is decomposed into two parts, 
the matter (T^'') and the electromagnetic (T^'^j) parts: 



^ — -'-m + -'-EM- 



(8) 



The electromagnetic part is given by 

nil = ^ [f^F^" - \9'''F,,FP'^^ . (9) 

In general the stress-energy tensor for the accretion disk 
is determined by mass density, internal energy, pressure, 
viscosity, radiative transfer and etc.|l7j. Since we are 
interested in the accretion dominated by the Poynting 
flux, the two-dimensional accretion disk is assumed to be 
non- viscous, cool, and non-radiative such that the matter 
part is given by 



(10) 



where pm is the rest-mass density and u'^ is the four 
velocity of the accreting matter. 

From the conservation of stress-energy tensor, one can 
obtainQ 



{drUo)M+ + 2TrrK^{-aE^ + pmB^) = 0, 
{drU^)M+ + 2TTrK^mB^ = 0, 



M4 



+ u 

-g-^ — {drU^) 



n, 



drU^ 



0. 



(11) 
(12) 



(13) 



where the angular velocity of the disk is given hy fin = 
u"^ /vP and the mass accretion rate is given by 



(14) 



Ea. (|ll() and H12|l correspond to the energy and angu- 
lar momentum conservation respectively in the station- 
ary and axisymmetric setting in this work. One can see 
that the accretion is determined by the Poynting flux. 
The last equation is obtained from the radial component 
which essentially determines the orbital motion of the 
disk. Since = 0, the 0-component of the conserva- 
tion equation does not contribute to the accretion flow. 
In the absence of external fields, the angular velocity is 
determined as 



D 



drUQ 



(15) 



which is one of the characteristics of the Keplerian orbit. 
Hence deviations oi from the Keplerian one is natu- 
rally expected in the Poynting flux dominated accretion 
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disk. One of the purposes of this work is to see how it 
depends on the magnetic field. 
In the force- free hmitlla, 



we get 



E 



(16) 



(17) 



in the Schwarzschild background. $1^? is the angular ve- 
locity of the magnetic surface which rotates rigidly in an 
axisymetric and stationary state . Then the accretion 
equations driven by the Poynting flux can be rewritten 
as 



{drUo)M+ + rvDVLpB^B^ = 0, 



M4 



27rr2- 



drUo 



-nl -11=0. 



(18) 
(19) 



(20) 



From ea. H18() and (|19|l . one obtains an interesting relation 



drUo 
drU^ 



(21) 



which has no explicit dependence on the field configura- 
tion. It seems to imply that Jlj^ is determined essentially 
by the dynamics of the disk. However one should note 
that the dynamics itself is governed not only by the grav- 
ity but also by the electromagnetic field as well. 



B. Stream equation 

The configuration of the ordered magnetic field around 
the accretion disk has been discussed both in analyti- 
cal and numerical studies. The developments of the or- 
dered magnetic field in the disk and the Poynting outflow 
from the disk have been studied numerically by many 



authors[T3IIllI3- 

For the strong enough electromag- 
netic field around the compact object, the force- free mag- 
netosphere can be established. In the non-relativistic for- 
mulation, Blandford*^ suggested an axisymmetric and 
stationary electromagnetic field configuration around an 
accretion disk on the equatorial plane. The poloidal 
field configuration for a black hole is known to satisfy a 
second order elliptical differential equation called Grad- 
Shafranov eauation|l9j or a stream eauation|2Clj| for the 
stream function \1/ and current /. The poloidal and 
toroidal components of the magnetic field can be writ- 
ten in terms of 5* and / respectively: 



B^ = 



2'K'UD 

21 



B'f' = - — . (22) 
zua 

Possible types of solutions in a force-free magneto- 
sphere has been discussed recently in the relativistic 
formulation]^ l2l| . 

In a fiat background, the stream equation becomes sim- 
pler one given by 



a^^- + ( -^de'f ] - np sin^ 619^ (r^l^F^r*) 



sin 9 



dl 



-flpsmede {sineflFde'i') = -167r^/— (23) 

BlandfordQ suggested two types of solution of ea. H23|) 
depending on the types of surface, where magnetic field 
lines rotating with angular velocity flp lie on. In a up- 
per hemisphere, the magnetic field of paraboloidal type 
cutting the disk on the equatorial plane is given by 



B' = 
B' = 
B^ = 



C 



2r{l + nlr^{l - cosOyy/^' 



-C(l - cos 6*) 



2rsin6'{l + nlr^{l - cos6')2}i/2' 
-Cflpil - cos 61) 



sin6l{l + 17|r2(l-cos6l)2}i/2' 
and the solution of hyperboloidal type is given by 



(24) 
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Boi?^(l - cos 61)2 
{sin^ 6'i?2 _ 7-2(1 - cos6l)4}{i?2 _ l]|r4(l - cos6l)4}i/2 

(1 - cos 6l){sin2 61 - (1 - cos 6')2}l/2 + {^2 _ ^2 (--^ _ 6l)2}l/2 COS 61 
Soi?2(l- COS 61)2 

{sin2 6'i?2 - r2(l - cos6l)4}{i?2 _ Vfipr^il - cos6l)4}i/2 

(1 - COS 6l){sin2 61 - (1 - cos6')2}i/2cos6l/sin6' - {R^ -r'^{l- cosOfy/'^ sint 

-BoRril - cos6')2l7^{sin2 OR^ -r^{l- cos 0)^} 
sin6l{sin6'2i?2 _ ^2(1 _ cos6')4}{i?2 _ nj,r^{l - cos6')4}i/2 ' 



(25) 



where C and represent the scale of field strength and 
R corresponds to the shape parameter for the hyper- 
boloidal configuration. It is interesting to note that Qp 
is an arbitrary function of r and 9. 

For the accretion flow, the magnetic field on the equa- 
torial plane, 9 = tt/2, is relevant. Among the accretion 
equations the implications of ea. (|2()|l are quite different 
for those two configurations. For the hyperboloidal type, 
B"" = on the equatorial plane in ea. H25ll and ea. H2U|) 
reads 



in the frame work of a simplified two-dimensional disk 
model on the equatorial plane. The physical variables 
for this simple system;30j are vP,u^,u'^, 5*, / and M+. 
Among those variables only M+ is independent of the 
disk location and we take it as an input parameter in 
this work. Then we can solve the equations for this sim- 
ple disk model to see the effect of the magnetic field on 
the accretion flow as well as to find the field configuration 
subject to the steady state accretion fiow driven by the 
magnetic field itself. 



0, 



(26) 



which may be considered to be the perfect-conductor 
boundary condition on the accretion disk. For the 
paraboloidal configuration, ea. (|20|l can be written as 



r^nl - 1) 



(27) 

Hence one cannot expect flu = flp unless vp{— rflp) = 
1 everywhere, which seems to be unphysical. It is natu- 
ral to suppose that the velocity of the magnetic field line 
vf{— rftp) does not exceed the speed of light. Then 
the r.h.s of the equation is positive. Since the term 
M+{vP lu'^)drU^ of the equation is negative, the term 
[VId — ^f] has to be negative such that 



^p > • 



(28) 



With a black hole at the center, however, the flat back- 
ground is no more valid and the field configuration is 
expected to be changed accordingly. Since there is no 
known analytic solution in the background of a black 
hole, we have to rely on the numerical calculation. The 
detailed numerical calculation will be discussed in the 
next section particularly for the paraboloidal configura- 
tion in the Schwarzschild background. 

In summary, we have a simple picture of accretion flow 
in an axisymmetric, stationary and force-free magneto- 
sphere, which is described by the above four equations 



III. NUMERICAL RESULTS FOR A 
PARABOLOID-TYPE CONFIGURATION IN 
SCHWARZSCHILD BACKGROUND 

In a fiat background, we can find / and 5" correspond- 
ing to the paraboloidal type configuration given by 



I{X) 



±cnFX 



2{i + filx^y/^' 



ttC 



where 



dX (H- 01x2)1/2 ' 



X = r(l ±cos6'). 



(29) 



(30) 



(31) 



The explicit forms for the magnetic fields are just those 
given in the previous section, ea. (|24|l . When fij? is con- 
stant, we can easily integrate d'^/dX of ea. (p?n|l to get 
explicit form of ^E* as 



*(X) = ^smh-\npX). 
ilp 



(32) 



In the Schwarzschild background, the stream equation 
is given by 



5 



a T 



r'' \sm9 J 

Q 1 fi™-2 T dl 



1 _ 2M 



(^l — ML 



(33) 



r 



Compared to the case with flat background we do not 
have analytical forms like ea. l29|l and (jSOJ. When flp 
and / are vanishing, ea. H33l) becomes 



dr 



1 - 



2M 



smf 



■df> 



= 0. (34) 



One of the solutions, ^fg, suggested by Blandford and 
Znajek [H is given by 



(35) 



where 



X EE r(l=Fcos6l) + 2M(l±cos6l){l-log(l±cos6l)}, (36) 

and upper/lower sign refers to the upper/lower hemi- 
sphere. One can easily see that it reduces to the 
paraboloidal field lines for M ^ or r ^ M. The up- 
per (lower) sign corresponds to the northern (southern) 
hemisphere. It is interesting to note that the discontinu- 
ity of on the disk surface (6* — 7t/2) defines the toroidal 
surface current K'^ 



C 



(37) 



which is the same as in the flat background^. 

For non- vanishing f2 p and J, we consider a solution for 
which the shape of the magnetic surface are the same as 



the magnetic surface defined by 22'|. Thus 5",/ and 
Hp are assumed to depend on X. We suppose that the 
derivative of 5* has the same form as in the flat back- 
ground, ea. H3U|) . That is, we take an Ansatz such that 



dx {i + nix^y/^' 



(38) 



On the equatorial plane, X = r + 2M, we can construct 
magnetic fields on the disk surface as follows: 



2r(l + f^|,X2)i/2' 
-C 



2r(l + f^|,X2)i/2 



1/2 



r I 



-2/ 



2M 
r 



-1/2 



Using ea. l|^ . eq.ljSSI) can be written as 



AttCVLp 



(1 + 17|,X2)3/2 (X - 2Mf 



where 



(39) 



(40) 



(Kip 

Ix 



{X ~ 3M)(4M^ - 6M'^X + 5MX^ - X^) + Qp {lOM^ - ISM^X + MIX^ - X"^ + 3M^X{X - 2Myn%}{Al) 



It can be considered to be an equation for / on the disk 
surface(0 = '7t/2) provided fJi^ is given. In contrast to the 
flat background case, we do not have the analytic form 
for / analogous to eq. . 

The accretion equations can then be written as follows. 



{drUo)M+ + 



{i + nlx^' 

C2 

{i + nix^ 



-1 = 0, 



■1 = 0, 



(42) 
(43) 



(44) 

where a dimensionless quantity / for the upper hemi- 
sphere is defined by 

i^ii + nWf^'^ . (45) 

Four basic equations, ea. H4U|) - ea. H44|) . governing the 
accretion flow under the influence of paraboloidal-type 
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configuration are solved numerically j^^. We get the ra- 
dial variations of Uq, m^, / and Qp on the disk with 

dimensionless variable M+ — M^/C'^ as an input pa- 
rameter in this work. It should be noted that the di- 
mensionless accretion rate parameter M-j- is the relative 
magnitude of the accretion rate to which determines 
the strength-squared of the magnetic field. 

For numerical calculations, we need a set of initial val- 
ues to start with. We take an initial point tq far from 
the horizon such that we can make use of the Newtonian 
approach in the initial step of the numerical calculation. 
We take to be much greater than inner most stable 
radius 

ro > 6M. (46) 

In this asymptotic region we can take initial values as 
those for the flat background approximation. At first, 
we choose initial values for m° and u'^ to be 

u' = ul, u'^ = vu%, (47) 

where and are those for the Keplerian orbit. We 
take the initial value of u'^ to be smaller than that of the 
Keplerian orbit to prevent from becoming an imagi- 
nary number, rj < 1. Plugging the initial values of m'^ and 
u'^ into ea. (|43l) and eg . (|44|l . the initial value offlp can be 
determined. And the initial value of / is then determined 
by making use of of flat-background relation in eq. H29() : 

J = ^npr. (48) 

After fixing these initial values, it is straightforward to 
execute the numerical calculation. 

For numerical calculation, we consider two cases of the 

accretion-rate parameters, [a] M+ = 0.7 and [b] = 3 
and take initial point to be ro = 50Af. To get some 

physical idea on the values of Af+, let us take C to be 
2.78x10"'^^ gauss-cm for a solar mass black hole as an 
example. Then the magnetic field strength becomes ~ 
10^^ gauss at r = ro which might be be compatible with 
the astrophysical phenomena accompanied with strong 
magnetic field, like magnetar and GRB 0|. Then the 
accretion rates are given by A/+ = 1.8 x lO^^g/s and 

7.7 X lO^^g/s for M+ = 0.7 and 3 respectively. 

We continue the numerical calculation inwardly up to 
r = 6M which is the radius of the innermost stable cir- 
cular orbit in the Schwarzschild black hole at the center. 
As shown in Figure ^ the angular velocities of the mag- 
netic field lines(rii?) are found to be different from either 
the Keplerian angular vclocity(r2/f ) or the disk angular 
velocity(ri£)). It implies that the two-dimensional disk 
model may not be a good approximation for the accre- 
tion disk of perfect conductor. As expected in the pre- 
vious section for the velocity of the magnetic field line 
VFi= r^lp) less than the speed of light, Hp is found to 
be larger than flpi. 



For a smaller accretion rate, M+ = 0.7, flo changes its 
sign at r ~ 20 AI as shown in Figure ^a]. It corresponds 
to changing its sense of rotation. It is because for a given 
strength of the magnetic field, the slope of angular mo- 
mentum {drU,p) is proportional to the inverse of Af_|_ as 
can be seen in ea. ()19f) . Hence the angular momentum 
of the accretion fiow decreases more rapidly for smaller 
Af-|_ and there is a chance of reverting the direction of 

rotation. But for larger A/+, the slope of change is not 
high enough to change the sense of rotation. Figure ^b]. 

We obtain also the numerical solutions for and / of 
the stream equation to get the electromagnetic field on 
the disk. It is found that the strength of the magnetic 
field is increasing as r goes near the inner edge as shown 
in Figure [3 One can see that the axial component B'^ 
is sensitive on while and B^ are not sensitive to 

M -|_ . The surface currents are defined as discontinu- 
ities of the magnetic fields across the disk and they show 
similar behaviors as those of the magnetic field. Since 
the magnetic field as well as the surface current are in- 
creasing as r gets smaller, it is naturally expected that 
the Poynting flux increases substantially as r approaches 
to the center as shown in Figure 13 Although the nu- 
merical calculation does not go beyond r < 6A/, it may 
indicate a possible electromagnetic jet structure near the 
inner edge of the disk. 



IV. DISCUSSION 

In this work, we investigate the accretion flow domi- 
nated by Poynting flux using a simplified model of a two- 
dimensional disk. The two-dimensional accretion disk 
in a force-free magnetosphere is described by the basic 
four equations, three accretion equations and one stream 
equation. We observe that the condition of the stationary 
accretion flow (constant M+) gives rise to nontrivial con- 
straints on the form oiQp and on the configuration of the 
ambient electromagnetic field. It is found that it depends 
strongly on the types of field. In the fiat background 
limit, we find that the hyperbloidal type configuration 
suggested by Blandford can accommodate the solutions 
with flp = Qpi, which is consistent with an accretion 
disk of perfect conductor. However for the paraboloidal 
type configuration, Oi? is expected to be greater than 
flo- Hence one can guess that a two-dimensional flow 
with paraboloidal type configuration might be relevant 
for an accretion disk not in a perfect-conducting state. 

To get more details on the dynamics of the accre- 
tion flow, we try to find numerical solutions particularly 
for the paraboloidal type configuration in Schwarzschild 
background. The overall dynamics of the disk is found 
to depend essentially on the ratio of the accretion rate to 

the field strength-squared, Af+. When it is sufficiently 
small, the inner part of the accretion disk is found to 
be changing its sense of rotation to match the angular 
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momentum balance with respect to the Poynting flux. 

One of the interesting results of this work is to demon- 
strate clearly how fij? is different from either flo or 
under the influence of the strong magnetic field. 
Since the numerical calculation in this work is extended 
only up to the inner stable radius for the Schwarzschild 
black hole, the relativistic effect is found not to be dom- 
inant such that the discussion in a flat background limit 
does not change substantially. Hence the observation of 
— for hyperboloidal type configuration in a two- 
dimensional accretion fiow is expected to be valid even 
with Schwarzschild black hole at the center. 

We demonstrate in this work that how the accretion 
disk provides relevant boundary conditions, which is im- 
portant to find the solutions of the stream equation. 
In this work it is possible basically because we have a 
mathematically well-defined disk surfaces in this simpli- 
fied model as in the case of horizon. However it would 
be interesting to see whether any realistic accretion disk 



model can provide the relevant boundary conditions for 
the stream equation, to which our analysis can be ap- 
plied. 

In order to simulate the strong relativistic effect it 
is necessary to extend the present study of the two- 
dimensional accretion fiow to that in the Kerr back- 
ground. However as the inner edge of the accretion disk 
gets very near/or beyond the ergo-sphere for a rapidly 
rotating black hole, the magnetic c oupling between the 
black hole and the accretion disk |23.l25ll26| becomes im- 
portant and the simple picture in this work should be 
modified substantially. 
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FIG. 1: The angular velocities of the magnetic surface Qf, accretion disk flo and Keplerian orbit VIk are shown in solid, dotted 
and dashed lines respectively for two cases of accretion rate: [a] M+ = 0.7 and [b] M+ = 3. 
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FIG. 2: The strength of the magnetic field in the upper hemisphere: the axial component B'^, and B^ are represented by 
the dotted, solid and dashed lines respectively. 



